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A  conical  surface  roughness  model  applicable  to  particle  simulations  has  been  developed.  The 
model  has  been  experimentally  validated  for  channel  flows  using  helium  and  nitrogen  gases  at 
Reynolds  numbers  from  0.01  to  10  based  on  inlet  conditions.  To  efficiently  simulate  gas-surface 
interaction,  molecular  collisions  with  the  actual  rough  surface  are  simulated  by  collisions  with  a 
randomly  positioned  conical  hole  having  a  fixed  opening  angle.  This  model  requires  only  one 
surface  parameter,  average  surface  roughness  angle.  This  model  has  also  been  linked  to  the 
Cercignani-Lampis  scattering  kernel  as  a  required  reference  for  use  in  deterministic  kinetic 
solvers.  Experiments  were  conducted  on  transitional  flows  through  a  150pm  tall,  1cm  wide, 

1.5cm  long  microchannel  where  the  mean  free  path  is  on  the  order  of  the  roughness  size.  The 
channel  walls  were  made  of  silicon  with:  (i)  polished  smooth  surfaces,  (ii)  regular  triangular 
roughness,  and  (iii)  regular  square  roughness  with  characteristic  roughness  scales  of  <1  pm,  11pm, 
and  29pm  respectively.  For  the  triangular  roughness,  mass  flow  reductions  ranged  from  6%  at  the 
higher  stagnation  pressures  tested  to  25%  at  the  lower  stagnation  pressures  tested  when  compared 
to  the  smooth  channel. 
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I. 


INTRODUCTION 


Surface  structure  and  roughness  plays  a  role,  to  some  extent,  in  all  flow-surface 
interactions  regardless  of  speed  or  scale.  Traditionally,  surface  roughness  is  considered 
important  in  low  speed  micro-flows,  where  the  surface  to  volume  ratio  is  large  and  boundary 
conditions  dominate  the  flow  [1].  Since  the  first  report  was  published,  almost  a  century  ago  [2], 
a  substantial  amount  of  research  on  the  impact  of  surface  roughness  on  gas  flows  in  ducts  has 
been  accomplished.  For  instance,  the  mass  flow  degradation  due  to  surface  roughness  was 
shown  numerically  and  experimentally  by  Davis  et  al.  [3].  More  recently,  significant  effect  of 
surface  roughness  on  pressure  and  heat  transfer  in  circular  tubes  and  channels  was  demonstrated 
experimentally  in  Mala  and  Li  [4].  In  the  experiments  by  Turner  et  al.  [5]  however,  the  surface 
roughness  was  found  to  have  a  small  effect  on  gas  pressures  inside  a  channel  for  several  pressure 
ratios.  In  addition  to  low  speed  flows,  surface  roughness  has  been  seen  to  cause  a  noticeable 
impact  on  some  high-speed  flows,  such  as  micro-nozzle  expansions  [6]. 

The  applicability  of  continuum  numerical  approaches,  such  as  those  based  on  the  solution 
of  Navier-Stokes  equations  [7],  is  limited  for  flows  where  the  roughness  size  is  smaller  than  the 
gas  mean  free  path.  There  has  been  work  recently  on  slip  conditions  to  which  overcome  this 
limitation  [8]  In  most  of  the  aforementioned  applications,  the  characteristic  roughness  size  is 
smaller  than  1pm,  and  is  usually  on  the  order  of  or  smaller  than  the  gas  mean  free  path.  Being 
essentially  a  kinetic  effect,  the  interaction  of  gas  molecules  with  such  small  surface  roughness 
requires  a  kinetic  numerical  approach  in  order  to  properly  address  the  gas-surface  interaction  at 
the  level  of  molecular  velocity  distribution  functions.  Kinetic  modeling  of  surface  roughness 


Approved  for  public  release;  distribution  is  unlimited 


2 


effects  [9]  is  hindered  by  a  number  of  problems,  from  the  accurate  representation  of  complex 
shapes  and  roughness  patterns  to  difficulties  of  experimental  validation. 

The  fonner  has  been  approached  by  a  number  of  models,  different  in  their  numerical 
complexity  and  physical  accuracy.  Among  simplified  models  are  the  Cercignani-Lampis  (CL) 
kernel  [  10][  1 1  ]  and  cone  model  [  12][  13].  In  the  CL  model,  rough  surfaces  may  be  simulated  by 
specifying  the  tangential  accommodation  coefficient  at  larger  than  unity  [14].  A  similar  idea 
was  applied  by  Ketsdever  et  al.  [15]  where  an  anti-specular  reflection  was  introduced.  In  that 
model,  a  fixed  fraction  of  molecules,  a,  is  reflected  so  that  the  reflected  normal  and  tangential 
momentum  is  the  same  magnitude  as  the  incoming,  but  with  opposite  directions.  The  models 
[10]  and  [15]  are  easy  to  implement  and  may  be  used  with  particle  as  well  as  deterministic 
kinetic  solvers,  but  they  do  not  have  a  clear  correlation  between  the  parameters  of  the  models 
and  the  physical  properties  of  the  rough  surface.  This  drawback  was  overcome  for  particle 
simulations  in  the  cone  models  [  12] [  1 3]  where  the  surface  is  either  represented  by  a  number  of 
virtual  cones  or  pre-computed  distribution  functions  are  used  at  the  processing  stage.  These 
models  refer  to  the  actual  average  angle  of  the  rough  surface,  although  it  is  hard  to  use  a 
prescribed  angular  distribution,  and  a  special  cut-off  has  to  be  used  limiting  the  number  of  cones. 
Another  approach  is  to  explicitly  specify  the  rough  surface  as  a  set  of  randomized  or  regular 
structures  [16].  While  this  is  a  reasonable  approach  for  simple  geometries,  it  may  be  difficult  to 
use  for  complex,  especially  three-dimensional,  surfaces  and  roughness  structure.  Among  the 
more  sophisticated,  although  difficult  to  implement,  roughness  models  is  the  fractal  model 
[  17][  18]  and  an  analytic  model  [19]. 

The  experimental  validation  of  rough  surface  models  is  hampered  by  requirements  of 
high  precision  and  small  scale,  both  of  which  increase  the  difficulty  of  experimental 
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measurements,  as  well  as  the  superposition  of  different  physical  effects  [20],  Examples  of  related 
work  on  surface  roughness  are  [14]  where  flow  rates  through  short  channels  were  measured  and 
calculated  with  a  Monte  Carlo  method,  and  [21],  where  triangular  roughness  in  channel  flows 
was  studied  numerically  and  experimentally  for  large  Knudsen  numbers.  Adding  to  the 
difficulty  of  numerical  validation  is  the  complexity  of  relating  the  measured  macro  parameters  to 
the  micro-phenomenon  of  interest  to  the  investigation.  Although  changes  in  macro-parameters, 
such  as  pressure  distribution  and  mass  flow,  can  lead  to  changes  in  micro-parameters,  such  as 
accommodation  coefficient,  these  macro-parameters  do  not  necessarily  indicate  specific  reasons 
for  the  change  in  accommodation.  Changes  in  angle  of  incidence  [22],  gas  temperature,  surface 
temperature,  surface  material,  and  surface  cleanliness  can  all  lead  to  changes  in  accommodation. 
To  complicate  matters,  these  parameters  can  be  difficult  to  quantify  experimentally.  In  addition 
to  the  change  in  accommodation  due  to  the  surface,  a  long  microchannel  flow  may  be  affected  by 
variation  of  accommodation  with  local  gas  pressure.  In  most  cases,  the  probability  of  surface 
sticking  /  accommodation  decreases  as  the  availability  of  surface  sites  decrease.  Lundstrom  [23] 
observed  the  surprising  result  that  Knudsen  flow  diffusion  increased  with  backing  pressure.  This 
was  hypothesized  to  be  due  to  the  variation  in  sticking  with  local  gas  density  and  the  tendency  of 
molecules  to  scatter  more  specularly  off  an  absorbed  wall  gas  then  off  the  surface. 

This  study  consists  of  a  numerical  investigation  and  experimental  validation  of  the  effects 
of  surface  roughness.  Specifically,  the  scope  will  be  limited  to  two-dimensional  helium  and 
nitrogen  microchannel  flows  with  well-characterized  surface  roughness.  These  flows  will  be  in 
the  low  Reynolds  number  regime,  based  on  entrance  conditions,  with  surface  rarefaction  from 
near  free  molecular  to  near  continuum,  based  on  surface  roughness.  A  simple  model  of  the 
surface  roughness  that  maintains  the  detailed  balance  and  is  applicable  to  direct  simulation 
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Monte  Carlo  (DSMC)  method  has  been  developed,  and  its  connection  to  the  Cercignani-Lampis 
scattering  kernel  is  established.  This  new  model  is  more  efficient  in  both  setup  cost  and 
simulation  time  than  the  direct  representation  of  the  physical  geometry.  It  does  not  exhibit  the 
shortcomings  seen  in  previous  simple  surface  roughness  models  while  being  versatile  enough  to 
incorporate  either  an  average  surface  angle  or  a  distribution  of  surface  angles,  as  well  as  arbitrary 
roughness  coverage. 


II.  NUMERICAL  APPROACH 

A.  Conical  Model 

Molecular  collisions  with  a  rough  surface  can  be  separated  into  two  stages:  (1)  a 
molecule  strikes  the  wall  at  a  particular  point  and  (2)  the  molecule  experiences  one  or  more 
collisions  with  the  surface  roughness  feature  before  leaving  the  surface.  When  the  characteristic 
roughness  size  is  significantly  less  then  the  mean  free  path  of  the  gas,  no  intermolecular 
collisions  occur  during  the  second  stage.  The  time  taken  for  this  stage  is  much  smaller  then  the 
mean  collision  time,  and  the  distance  traveled  is  on  the  order  of  the  roughness  size.  Therefore, 
for  numerical  simulation,  the  time  and  distance  can  be  ignored  and  it  is  sufficient  to  specify  only 
the  reflected  velocity  of  the  molecule  after  the  collision. 

To  represent  this  two-stage  process  in  an  efficient,  yet  representative  way,  a  conical 
model  is  used  to  select  the  reflected  velocity  of  the  collided  molecule.  At  the  point  of  collision, 
the  rough  surface  is  represented  as  a  randomly  positioned  virtual  conical  hole  with  fixed  opening 
angle  />'  and  height  h.  The  molecule  is  assumed  to  enter  the  cone  through  the  base  at  point  A,  as 
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seen  in  Figure  1.  The  collision  point  for  the  molecule  and  the  cone  wall  is  then  calculated  at 
point  B.  Here  the  molecule’s  reflected  velocity  is  selected  assuming  a  diffuse  reflection, 
although  a  different  reflection  model  can  be  used.  The  process  of  calculating  subsequent 
collision  points  and  determining  post  collision  velocities  is  continued  until  the  molecule  leaves 
the  cone.  The  point  of  origin  for  the  reflected  molecule  is  assumed  to  be  the  same  as  the 
incoming  point  of  collision;  all  movement  in  the  virtual  cone  is  only  a  means  for  determining  the 
reflected  velocity  of  the  colliding  molecule.  Therefore,  the  actual  value  of  h  is  not  important 
(e.g.  h=  lm  can  be  used).  In  this  algorithm,  there  are  no  uncertainties  associated  with  cutoff  for 
very  long  traveling  distances  such  as  those  in  [12].  Due  to  surface  shape  simplification,  this 
algorithm  requires  only  one  parameter,  the  average  surface  slope  (represented  as  the  cone  angle 
(3),  making  it  easy  to  implement  in  DSMC. 

To  study  the  influence  of  a  rough  wall,  represented  by  the  conical  model,  on  flow 
properties,  examine  how  the  velocity  distribution  of  incident  molecules  is  transformed  by 
collisions  with  the  wall.  Consider  a  distribution  of  incoming  molecules  from  an  equilibrium 
flow  characterized  by  a  Mach  number  M  and  directed  along  the  x-axis,  assuming  the  rough 
surface  to  be  in  the  xy  plane.  This  scenario  has  been  simulated  using  the  algorithm  described 
above  and  the  distribution  of  the  angle,  (p ,  between  the  reflected  tangential  molecular  velocity 
and  the  flow  direction  is  plotted  in  Figure  2.  Note  that  the  maximum  of  the  distribution 
corresponds  to  the  direction  opposite  to  the  flow,  which  is  generally  expected  to  increase  flow 
resistance. 

Verification  for  the  conical  surface  model  has  been  performed  for  a  two-dimensional 
thennal  bath  to  assure  adherence  to  detailed  balance.  The  thermal  bath  constituted  helium  gas 
initially  heated  to  1000K  within  a  2D  box  with  a  surface  temperature  of  300K.  Gas  surface 
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collisions  were  simulated  by  the  conical  roughness  model  with  a  cone  angle  of  [3=45  deg.  The 
temperature  relaxation  time  inside  the  test  box  is  illustrated  in  Figure  3  for  a  cross  section  along 
the  centerline.  As  particles  collide  with  the  surface,  the  temperature  decreases  from  its  initial 
value  to  the  equilibrium  value  of  300K.  By  coming  to  a  stable  equilibrium  the  simulation  proves 
adherence  to  the  detailed  principle. 

B.  CL  Scattering  Kernel 

The  above  stated  cone  model  technique  is  valuable  for  particle  simulations,  but  does  not 
have  a  direct  applicability  in  deterministic  kinetic  solvers  since  it  is  based  on  individual  particle 
velocities  and  not  on  distribution  functions.  It  is  reasonable  to  establish  a  connection  between 
the  conical  model  and  a  kernel  that  is  directly  applicable  to  a  kinetic  solver.  One  such  kernel,  the 
Cercignani-Lampis  (CL)  [10], [24]  scattering  kernel,  can  be  used  in  kinetic  solvers  but  lacks 
correlation  between  the  physical  surface  and  the  main  input  parameter  for  the  kernel.  The 
conical  surface  model  can  be  used  to  relate  this  parameter,  the  tangential  momentum 
accommodation  coefficient  ah  to  an  average  surface  roughness  angle. 

To  facilitate  this  relation  it  is  necessary  to  create  an  analytic  expression  for  the 
distribution  of  (p  as  a  function  of  at.  Assuming  the  temperature  of  the  gas  and  the  wall  are  equal 
and  velocities  are  normalized  by  yjm/ kT ,  the  distribution  functions  of  x  and  y  velocity 
components  of  incident  molecules  are 


respectively.  After  the  CL  transfonnation,  the  x  component  of  the  reflected  velocity  is  equal  to 
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(3) 


Vrefl  =  1  a'  vinc+vCL 


VF 


a. 


where  the  distribution  function  of  vx  is  given  by 


/f(v)  = 


2;r(l-|l-a,|) 


exp 


2;r(l-  1  -at ) 


(4) 


Here,  at  is  the  CL  parameter  that  corresponds  to  the  tangential  accommodation  coefficient,  0<  at 
<2.  a,  =1  corresponds  to  zero  average  tangential  momentum  relative  to  the  surface;  the  larger  the 
deviation  between  a,  and  unity,  the  larger  the  average  momentum  of  the  reflected  molecules. 
Note  that  Vf  is  the  sum  of  the  two  independent  normally  distributed  components,  therefore,  its 


distribution  function  can  be  written  as 


/(v7)  =  — exp 
2  n 


,2V 


refl 


~ M ^ 


]-a, 

VI1  ~a> 


(5) 


The  distribution  function  of  the  y  velocity  component  does  not  change  during  the  CL 
transfonnation.  Note  that  tan(^)  =  |v'',/7  |/v'"/z  ,  so  for  0  <(/)  <  k/2  the  distribution  function  of  (j) 


can  be  obtained  with 


ou 

/(tan  <f))  =  J ujrxefl (ux)fryefl (ux  tan  </>)dux 


(6) 


and  for  n/ 2  <(j)<7i 


/(tan  =  uxfxefl  (ux)fryefl  (ux  tan  (f))dux 


(7) 


Finally,  the  reflected  velocity  distribution  function  can  be  written 

L  (</>)  =  2/(tan  (/)){ 1  +  tan2  tj>) 


(8) 
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_f2 

ft (</>)  =  —+  —7=^ exp  (- 1 2  sin2  </>)(!  + erf  (t  cos </>))  (9) 

7t  v  n 

where  t  =  M y[y  {l-a^j ^|l-a,| 

Parameter  at  of  the  CL  transfonnation  can  be  found  by  least  square  fitting  of  (f> 
distributions  obtained  using  the  cone  reflection  algorithm  described  in  section  II: A  as  seen  in 
Figure  2.  Figure  4  shows  the  value  of  at,  found  by  a  least  square  fit,  for  M=0.3  and  different 
cone  angles  [3.  As  expected,  at  is  maximum  at  some  intermediate  value  of  the  cone  angle.  Since 

at  is  larger  than  1,  the  average  tangential  momentum  of  the  reflected  molecules  points  in  the 
opposite  direction  with  respect  to  the  average  tangential  momentum  of  the  incoming  molecules. 

The  value  of  at  only  weakly  depends  on  Mach  number,  as  seen  in  Figure  5,  where  at  is 
shown  as  a  function  of  Mach  number  for  two  different  opening  angles,  P  =120°  and  P  =66°.  This 
fact  facilitates  the  use  of  the  CL  model  in  simulations,  since  the  value  of  at  can  be  selected 
universally  depending  on  the  degree  of  roughness  of  the  surface,  and  not  on  the  flow  properties. 

C.  Simulation  Conditions 

In  order  to  validate  the  conical  roughness  model  and  its  link  to  the  CL  scattering  kernel, 
both  were  individually  incorporated  into  a  DSMC  solver  and  simulations  were  run  for 
comparison  against  experiment.  To  maximize  the  measurable  effect  of  the  different  surface 
models  on  the  flow,  the  simulations  were  conducted  on  a  long,  narrow  microchannel  with  an 
L/H  =100  and  an  assumed  infinite  width.  The  2D  capability  of  SMILE  [25]  was  used  in  all 
DSMC  computations.  This  system  utilizes  the  majorant  frequency  scheme  to  calculate 
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intermolecular  interactions  [26],  while  using  the  variable  hard  sphere  model  for  intermolecular 
potentials.  Four  surfaces  were  used:  (i)  a  fully  diffuse  smooth  surface  with  complete  energy 
and  momentum  accommodation,  (ii)  a  smooth  surface  utilizing  the  conical  surface  roughness 
model  with  an  opening  angle  of  66°  which  corresponded  to  an  earlier  experimental  setup 
concept,  (iii)  CL  model  with  a  tangential  momentum  accommodation  coefficient  of  1.045 
attained  using  the  conical  model  with  a  66°  cone  opening  angle,  and  (iv)  a  diffuse  surface 
consisting  of  about  900  triangles  representative  of  an  actual  experimental  shape. 

The  channel  length  and  half  height,  from  plane  of  symmetry  to  the  tops  of  the  surface 
features,  were  1.5  cm  and  75  pm,  respectively.  Helium  was  used  as  the  test  gas  with  stagnation 
pressure  from  200  Pa  to  8,000  Pa  for  all  four  surfaces.  The  Knudsen  number  is  defined  as 


1 

V2 n7id2L 


kT 

4lPnd2L 


(10) 


where  P  is  the  stagnation  pressure,  k  is  Boltzmann’s  constant,  T  is  the  stagnation  temperature 
300K,  d  is  the  molecular  diameter,  and  L  is  the  feature  size  (about  20  pm).  The  Knudsen  number 
for  these  tests  ranged  from  0.1  to  5.  Nitrogen  was  used  as  the  test  gas  for  only  surfaces  (i)  and 
(ii),  with  stagnation  pressures  from  200  Pa  to  6,000  Pa.  Again,  based  on  the  stagnation 
conditions  and  the  surface  roughness  size,  the  Knudsen  number  ranged  from  0.05  to  2.  A 
convergence  study  was  performed  to  assess  the  impact  of  the  inflow  boundary  conditions  on  the 
results  of  the  computations.  These  computations  were  perfonned  using  helium  at  6,000  Pa  and 
two  inflow  boundary  locations:  the  first  extending  1700  pm  from  the  channel  entrance  and 
1700  pm  from  the  symmetry  plane,  the  other  850  and  850  pm,  respectively.  A  Maxwellian 
distribution  function  with  zero  flow  velocity  was  assumed  at  these  boundaries.  No  impact  from 
the  boundary  location  was  found.  The  larger  domain  was  used  in  all  subsequent  computations, 
with  the  number  of  molecules  ranging  from  about  1.6  million  for  the  lower  pressures  to  about  10 
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million  for  the  larger  ones.  A  schematic  for  the  computational  domain  can  be  seen  in  Figure  6, 
where  the  outflow  boundary  condition  is  modeled  as  a  vacuum. 


III.  EXPERIMENTAL  APPROACH 

The  experiments  were  conducted  using  atomically  smooth,  triangular  patterned,  and 
square  patterned  silicon  textures.  Standard  microelectromechanical  systems  (MEMS)  processes 
were  performed  to  create  1  cm  by  1.5  cm  textured  areas  centered  on  2  cm  by  1.5  cm  chips.  The 
feature  width  for  both  square  and  triangular  textures  is  approximately  1 5  pm.  The  triangles  are 
approximately  1 1  pm  deep  with  an  etched  angle  of  54.7°  from  surface  parallel,  which  would  give 
an  opening  angle  of  approximately  70°,  very  close  to  the  numerical  simulation.  The  squares  are 
approximately  30  pm  deep,  both  geometries  can  be  seen  in  Figure  7.  Teflon  shims  were  placed 
between  the  chips  to  create  a  channel  1  cm  wide,  150  pm  high,  and  1.5  cm  long  such  that  gas 
flow  along  the  length  of  the  channel  was  perpendicular  to  the  textures  etched  in  the  silicon  chip. 
Along  this  length  there  were  approximately  900  features.  This  channel  was  secured  in  an 
aluminum  holder  and  dimensions  of  the  channel  were  measured  with  a  Scanning  Electron 
Microscope  before  and  after  testing  as  seen  in  Figure  8.  The  ratio  of  the  measured  open  area  of 
the  channel  to  the  open  area  of  the  numerical  simulation  was  used  to  linearly  correct  the 
experimental  mass  flow  data.  Typical  deviation  for  the  channel  height  was  3-5%.  Deviation  for 
the  width  was  less  then  1%. 

The  aluminum-Teflon-silicon  assembly  was  placed  in  the  wall  of  a  3550  cm3  plenum, 
which  acted  as  a  stagnation  chamber.  The  plenum  was  machined  with  two  gas  inlet  ports  and 
three  ports  for  pressure  measurements.  The  setup  was  operated  in  the  CHAFF  IV  vacuum 
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chamber  at  the  University  of  Southern  California,  which  is  capable  of  an  ultimate  pressure  of  10" 
6  Torr  with  working  pressures  no  higher  then  10‘4  Torr.  This  allowed  plenum  pressures  from 
approximately  1  to  100  Tom  with  ratios  of  stagnation  to  ambient  pressures  greater  than  105 
maintained  throughout.  Helium  or  nitrogen  test  gases  were  used.  Upstream  gas  flow  was 
measured  with  a  mass  flow  meter  and  regulated  through  a  needle  valve.  The  pressure  in  the 
plenum  was  monitored  with  several  differential  baratrons.  This  array  of  instrumentation  allowed 
measurements  of  pressure  and  mass  with  instrumental  precisions  below  1%. 

IV.  RESULTS  AND  DISCUSSION 

The  first  result  to  consider  is  that  of  the  numerical  simulation  using  the  smooth,  fully 
diffuse  surface  with  a  helium  test  gas.  Figure  9  shows  the  local  pressure,  normalized  by  the 
stagnation  value,  and  Mach  number  along  the  centerline  of  the  channel  for  plenum  pressures  of 
200  Pa  and  2000  Pa.  The  Knudsen  number,  based  on  channel  height  of  150  pm  instead  of  the 
feature  size,  for  the  200  Pa  case  increases  from  below  1  to  about  100.  As  seen  in  Figure  9,  the 
pressure  in  the  free  molecular  flow  drops  linearly  from  Po  =  200  Pa  to  0.0 lPo  through  the 
channel.  The  Mach  number  stays  below  0.05  for  the  majority  of  the  channel,  increasing  to  1  at 
the  exit;  the  Mach  number  for  the  Po  =  2000  Pa  case  is  higher,  as  expected.  This  allows  for  the 
use  of  a  constant  at  value  in  the  CL  model  along  the  length  of  the  channel. 

Figure  10  shows  the  ratio  of  local  pressure  for  the  conical  surface  simulation  divided  by 
the  local  pressure  for  the  smooth  surface  as  a  function  of  longitudinal  distance  down  the  channel. 
As  expected,  the  pressure  along  the  channel  increases  for  the  rough  surface  compared  to  the 
smooth,  since  the  rough  surface  impedes  the  flow.  This  effect  decreases  as  the  flow  becomes 
free  molecular,  at  which  point  the  pressure  ratio  down  the  channel  asymptotes  to  1  far  from  the 
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exit.  The  ratio  drops  below  1  at  the  point  where  the  local  mean  free  path  exceeds  the  channel 
height  allowing  molecules  to  “see”  the  exit  without  intermolecular  collisions.  Downstream  from 
this  point,  any  molecules  traveling  upstream  are  more  likely  to  leave  the  channel  after  colliding 
with  the  rough  surface  than  they  would  be  after  colliding  with  a  smooth  surface.  Upstream  from 
this  point,  molecular  collisions  are  likely  to  prevent  this.  Note  that  a  relatively  small  difference 
between  pressures  for  the  rough  and  smooth  surfaces  is  observed  for  much  of  the  channel  length, 
consistent  with  the  measurements  of  Turner,  et  al.  [5]. 

Calculated  helium  mass  flow  as  a  function  of  stagnation  pressure  is  listed  in  Table  1  for 
all  four  surfaces.  A  reduction  in  mass  flow  of  25%  at  the  lower  pressures  and  6%  at  higher 
pressures  is  seen  between  the  conical  surface  model  and  the  smooth,  fully  diffuse  case.  The  CL 
model  with  an  at  =  1 .045  matches  the  conical  model  to  7%  at  the  lower  pressures  and  less  then 
1%  at  the  higher  pressures.  The  conical  surface  matches  the  triangular  geometry  pattern  within 
3%  at  the  lower  pressures  and  less  then  1%  at  the  higher  pressures.  The  impact  of  the  surface 
roughness  is  larger  for  helium  then  for  nitrogen  for  any  given  pressure;  however,  when  plotted 
versus  Knudsen  number,  the  impact  is  similar.  Assumptions  made  for  the  cone  surface  model 
require  a  non-continuum  Knudsen  number  based  on  the  feature  size  of  the  roughness.  Although 
the  cone  model  agrees  favorably  with  the  direct  representation  over  the  full  range  tested,  it  is  not 
a  given  conclusion  that  the  cone  model  is  applicable  in  low  Knudsen  number  cases. 

There  is  good  agreement  between  the  experimental  smooth  silicon  surface  and  the 
smooth  fully  diffuse  simulation  as  seen  in  Figure  11.  The  experimental  uncertainty,  based  on 
standard  deviation,  varies  from  about  10%  at  Po  <  200Pa  to  less  then  1%  at  Po  >  10,000Pa. 
There  is  a  2%  under-prediction  by  the  DSMC  simulations  along  the  range  of  pressures  tested. 
One  of  three  effects  can  cause  this  deviation.  The  first  is  a  finite  specularity  that  may  have 
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existed  in  the  experiment  and  was  unaccounted  for  in  the  simulations.  The  second  is  an 
uncertainty  of  approximately  2%  in  the  area  of  the  channel  based  on  the  SEM  pictures.  The  third 
is  the  previously  mentioned  influence  of  surface  absorption  and  the  decrease  of  surface  sites  at 
higher  pressures,  also  creating  a  higher  physical  specularity  then  accounted  for  in  the  simulation. 
However,  relative  influence  of  these  effects  are  unknown.  It  is  important  to  note  that  the 
pressure  dependence  of  the  under-prediction  is  hard  to  determine  due  to  the  larger  experimental 
uncertainty  at  lower  pressures. 

There  is  also  good  agreement  between  the  experimental  triangular  surface  and  the  cone 
model  simulations  as  seen  in  Figure  12.  As  with  the  smooth  surface,  the  experimental 
uncertainty,  based  on  standard  deviation,  varies  from  about  10%  at  Po  <  200Pa  to  less  then  1%  at 
Po  >  10,000Pa.  At  the  higher  pressures,  there  is  less  then  2%  difference  between  the  experiment 
and  simulation;  however,  at  the  lower  pressures,  there  is  a  5%  over-prediction  by  the  DSMC. 
Again  both  the  uncertainty  in  the  experimental  channel  area  and  the  larger  statistical  scatter  at 
lower  pressures  could  contribute  to  this  difference. 

The  experimental  data  for  the  square  surface  features  can  be  seen  in  Figure  13  compared 
with  the  other  two  experimental  surfaces.  Molecular  interaction  with  the  vertical  walls  of  the 
square  surfaces  could  be  approximated  as  cones  with  opening  angles  (3=0°.  As  expected,  the 
square  data  falls  between  the  smooth  [3=180°  and  triangular  [3=70°  experimental  points. 
Although  the  square  surface  data  trends  towards  the  triangular  data  at  lower  pressures,  it  moves 
towards  the  smooth  data  at  higher  pressures.  This  could  be  attributed  to  the  influence  of  the  exit 
plane  on  the  flow  and  the  reduction  in  Knudsen  number  as  the  pressure  ratio  increases.  At  higher 
pressures  and  densities,  the  square  texture  acts  more  like  a  fully  diffuse  surface  then  a  series  of 
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roughness  features,  until  a  few  mean  free  paths  before  the  exit  plane,  ultimately  asymptoting  at 
the  smooth  surface  value. 


V.  CONCLUSION 

A  conical  surface  roughness  model  has  been  developed,  numerically  verified,  and 
experimentally  validated.  For  a  particle  simulation,  the  conical  model  reduces  both  setup  and 
simulation  time  over  direct  geometric  representation  without  sacrificing  accuracy.  For 
deterministic  solvers,  this  conical  model  gives  a  crucial  link  between  the  Cercignani-Lampis 
(CL)  scattering  kernel,  which  requires  a  defined  tangential  momentum  accommodation 
coefficient,  at,  and  a  physical  surface  characteristic,  the  average  surface  roughness  angle.  The 
conical  model  is  general  enough  to  be  used  for  surfaces  with  both  prescribed  roughness  features 
and  non-uniform  roughness  features.  In  the  latter  case,  a  mix  of  conical  opening  angles  and 
traditional  diffuse  interactions  can  be  used  based  on  the  surface  distribution. 

Numerical  simulations  using  the  conical  model,  the  CL  scattering  kernel,  and  a  direct 
geometric  representation  of  a  regular  triangular  surface  were  compared  against  a  smooth,  fully 
diffuse  model.  All  four  cases  were  simulated  with  the  DSMC  method  as  two-dimensional  flows 
of  helium  or  nitrogen  through  a  150pm  tall,  1.5cm  long  microchannel  with  Reynolds  numbers 
from  0.01  to  10  based  on  inlet  conditions.  The  direct  geometric  representation,  the  most  costly 
but  most  accurate  compared  to  experiment,  showed  a  decrease  in  mass  flow  between  25%  for  the 
lower  stagnation  pressures  simulated  and  6%  for  the  higher  stagnation  pressures  simulated.  The 
conical  model  matched  this  result  to  within  3%  at  the  lower  stagnation  pressures  and  less  than 
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1%  for  the  higher  stagnation  pressures  while  the  CL  model  matched  this  result  within  7%  and 
less  than  1%  respectively. 

Experiments  were  conducted  on  transitional  flows,  where  the  mean  free  path  is  on  the 
order  of  or  larger  than  the  roughness  size,  through  a  1cm  wide  channel  with  similar  height  and 
length  to  the  simulations  to  approximate  2-D  conditions.  The  channel  walls  were  made  of  silicon 
with:  (i)  polished  smooth  surfaces  (ii)  regular  triangular  roughness,  and  (iii)  regular  square 
roughness  with  characteristic  roughness  scales  of  <1  pm,  11  pm,  and  29pm  respectively. 
Experimental  results  validate  the  numerical  simulations,  which  match  to  within  2%  throughout 
the  range  tested  for  the  smooth  geometry  and  between  2%  and  5%  for  the  high  and  low 
stagnation  pressures  respectively  for  the  triangular  roughness  geometry.  The  experimental  mass 
flow  for  the  square  roughness  geometry  falls  between  the  smooth  and  triangular  cases  over  the 
entire  test  range  as  expected  based  on  p. 
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Pressure 

Pa 

Surface  (i) 
kg/s 

Surface  (ii) 
kg/s 

Surface  (iii) 
kg/s 

Surface  (iv) 
kg/s 

200 

4.132E-09 

3.188E-09 

3.440E-09 

3.093E-09 

650 

1.381E-08 

1.125E-08 

1.162E-08 

1.105E-08 

2000 

5.680E-08 

4.939E-08 

4.997E-08 

4.898E-08 

4000 

1.623E-07 

1 .484E-07 

1.502E-07 

1.483E-07 

6000 

3.094E-07 

2.899E-07 

2.893E-07 

2.896E-07 

8000 

5.098E-07 

4.810E-07 

4.830E-07 

4.810E-07 

Table  1 


Lilly 


Table  1  Numerical  results  for  mass  flow  as  a  function  of  stagnation  pressure  for  various  surface  models 
using  a  helium  test  gas. 
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Figure  1  Schematic  of  new  conical  model  reflection  algorithm. 


Figure  2  Distribution  function  of  the  angle  between  the  reflected  tangential  velocity  and  the  flow 
direction,  M=0.3,  (3=66°. 


Figure  3  Thermal  relaxation  of  helium  in  a  300K  box.  Time  is  in  mean  collision  times,  rx ,  at  1000K. 


Figure  4  Tangential  momentum  accommodation  coefficient,  at ,  as  a  function  of  cone  opening  angle  p, 
for  M=0.3. 


Figure  5  Tangential  momentum  accommodation  coefficient,  at  as  a  function  of  Mach  number  for  two 
cone  opening  angles. 


Figure  6  Schematic  of  the  computational  domain. 


Figure  7  Average  experimental  surface  roughness  feature  dimensions  for  a)  the  triangular  surface  and  b) 
the  square  surface. 
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Figure  8  SEM  pictures  of  a)  end  of  the  triangular  textures  (from  above),  and  b)  view  along  cleaved 
square  textures  (from  side). 


Figure  9  Normalized  pressure  and  Mach  number  as  a  function  of  distance  downstream  along  the 
centerline  for  two  stagnation  conditions. 


Figure  10  Pressure  ratio  as  a  function  of  distance  downstream  along  the  centerline,  conical  rough  surface 
to  smooth  surface. 


Figure  1 1  Comparison  of  the  experimental  smooth  texture  geometry  and  numerical  smooth  diffuse 
roughness  model  for  a)  helium  gas,  b)  nitrogen  gas. 


Figure  12  Comparison  of  the  experimental  triangular  texture  geometry  and  numerical  cone  roughness 
model  for  a)  helium  gas,  b)  nitrogen  gas. 


Figure  13  Experimental  results  for  the  textured  geometries  for  helium  gas. 
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